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OTI, 9UAMTY INSPECTED 



AERONAUTIC SYMBOLS 
1. FUNDAMENTAL AND DERIVED UNITS 


w 

9 

m 

I 


S 

S„ 

G 

b 

c 

A 

V 
2 
L 

V 
Do 
Dt 
Dp 
C 



Symbol 

Metric 

English 

Unit 

Abbrevia- 
tion 1 

Unit 

Abbrevia- 

tion 

Length 

l 

t 

F 

meter 

second 

m 

s 

kg 

foot (or mile) 

second (or hour) 

ft (or mi) 
sec (or hr) 
lb 

Force.. 

weight of 1 kilogram 

weight of 1 pound. 

Power 

P 

V 

horsepower (metric). 


horsepower 

a* 

fps 

Speed. 

/kilometers per hour 

(meters per second 

kph 

mps 

miles per hour 

feet per second 


2. GENERAL SYMBOLS 


Weight =mg 

Standard acceleration of gravity =9.80665 m/s 2 
or 32.1740 ft/sec 2 
W 

Mass=— 

Moment of inertia = to& 2 . (Indicate axis of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


v Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m -4 -s 2 at 15° O 
and 760 mm; or 0.002378 lb-ft~ 4 sec 2 
Specific weight of “standard” air, 1.2255 kg/m 3 or 
0.07651 lb/cu ft 


3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
Span 
Chord 

6* 

Aspect ratio, g- 
True air speed 
Dynamic pressure, gpl 72 


Lift, absolute coefficient C L —^ 

Drag, absolute coefficient G n =^ 

Profile drag, absolute coefficient Dd 0 =~^ 


D i 


Induced drag, absolute coefficient C D = ^ 


D„ 


Parasite drag, absolute coefficient Ddp—^ 
Cross-wind force, absolute coefficient Cc=^g 


G 


%n} 

it 

0 

B 


a 

e 

«o 

Oi a 

y 


Angle of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 


Reynolds number, p~~ where l is a linear dimen- 


sion (e.g., for an airfoil of 1.0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

Angle of attack 

Angle of downwash 

Anglo of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 

Flight-path angle 
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THE LAGRANGIAN MULTIPLIER METHOD OF FINDING UPPER AND LOWER LIMITS TO 

CRITICAL STRESSES OF CLAMPED PLATES 

By Bernard Bddiansky and Pai (’. Hu 


SUMMARY 

The theory of Lagrangian multipliers is applied to the prob- 
lem of finding both upper and lower limits to the true compressive 
buckling stress of a clamped rectangular plate. The upper and 
lower limits thus bracket the true stress, which cannot be exactly 
found by the differential-equation approach. The procedure for 
obtaining the upper limit, which is believed to be new, presents 
certain advantages over the classical Rayleigh-liitz method of 
finding upper limits. The theory of the lower-lim it procedure has 
been given by Trefftz but, in the present application, the method 
differs from that of Trefftz in a way that makes it inherently more 
quickly convergent.. It is expected that in other buckling prob- 
lems and in some, vibration problems the Lagrangian multiplier 
method of finding upper and lower limits may be advantageously 
applied to the calculation of buckling stresses and natural 
frequencies. 

INTRODUCTION 

Many important problems that cannot be exactly solved by 
the differential-equation approach and must therefore be 
analyzed by approximate methods arise in the buckling and 
vibrations of thin plates. The theory of Lagrangian multi- 
pliers can be a powerful tool in the analysis of many of these 
problems. The present paper presents the details of applica- 
tion as well as the fundamental principles of the Lagrangian 
multiplier method by demonstrating the use of the method 
to obtain both upper and lower limits to the true compressive 
buckling stress of a rectangular plate clamped along all edges. 

The procedure for obtaining the lower limit is similar to a 
method used by Trefftz (reference 1) and recently described 
by Reissner (reference 2). The present lower-limit method 
differs from that of Trefftz, however, in a way that makes it 
inherently more quickly convergent. The upper-limit pro- 
cedure, which does not appear to have been presented pre- 
viously, is simpler than the usual Rayleigh-Ritz method and 
may be expected to permit the computation of more accurate 
results with less labor. 

In a recent treatment of the problem of compressive buck- 
ling of clamped plates, extensive calculations of lower limits 
were made by Levy (reference 3) by means of a procedure 
equivalent to the Trefftz method. The results were esti- 
mated by Levy to be within 0.1 percent of the true results. 
In order to illustrate the methods of the present paper, upper 
and lower limits to the buckling stress of a square plate are 
computed to within 0.1 percent of each other; a positive 
check on the accuracy of Levy’s results is thus obtained. 
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SYMBOLS 

length of plate, in direction of stress 
width of plate, perpendicular to stress 
aspect ratio ( ajb ) 
thickness 
Poisson’s ratio 

Young’s modulus of elasticity 


Ef 


plate stiffness in bending ( 1 o( l—y 2 ) 

plate coordinate in direction of stress 
plate coordinate, perpendicular to direction of 
stress 

plate buckling deformation, normal to plane of 
the plate 

critical compressive stress, in ^-direction 

critical compressive stress coefficient in the 
formula <r x = k 

internal energy of deformation 
external work of applied stress 


F ourier coefficient of cos n -^ 


Fourier coefficient of cos rir ^ 
Fourier coefficient of cos 


a 

miry 


nicy 
cos , 
a b 


, j, m,n,p,q even integers 
', s odd integers 

Kronecker delta ( 1 if m — n : 0 if m^i 


mH |(>r ! re(i-)--km-8 -\ ( 1 -1 8 m \- S,J 

a, X, \j, y ( , y Lagrangian multipliers 

THEORETICAL BACKGROUND 

Rayleigh-Ritz method. — The Rayleigh-Ritz energy method 
for determining the critical stress of a thin plate consists of 
the following steps: 

(1) The deflection surface of the buckled plate is expressed 
in expanded form as the sum of an infinite set of functions 

1 
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having undetermined coefficients. In general, each term of 
the expansion must satisfy the geometrical boundary condi- 
tions of the problem. 

(2) The energy of the load-plate system is computed for 
this deflection surface and is then minimized with respect to 
the undetermined coefficients. 

(3) This minimizing procedure leads to a set of linear 
homogeneous equations in the undetermined coefficients. 
These equations have nonvanishing solutions only if the 
determinant of their coefficients vanishes. The vanishing of 
this stability determinant provides the equation that may he 
solved for the buckling stress. 

When the set of functions used is a complete set capable of 
representing the deflection, slope, and curvature of any possi- 
ble plate deformation, the solution obtained is, in principle, 
exact. Since, however, the exact stability determinant is 
usually infinite, a finite determinant yielding approximate 
results is used instead. 

Lagrangian multiplier method. — The Lagrangian mul- 
tiplier method follows the general procedure outlined for 
the Rayleigh-Ritz method, with but one outstanding change. 
The restriction in step (1) that the boundary conditions be 
satisfied by every term of the expansion is discarded and is 
replaced by the condition that the expansion as a whole sat- 
isfy the boundary conditions. This condition is mathemat- 
ically satisfied in step (2), during the minimization process, 
by the use of Lagrangian multipliers. 

The fundamental advantage of the Lagrangian multiplier 
method lies in the fact that, with the rejection of the neces- 
sity of the fulfillment of boundary conditions term by term, 
the choice of an expansion is much less restricted. In the 
elamped-plate compression problem, a simple Fourier expan- 
sion may be used instead of the complicated functions 
assumed in the Rayleigh-Ritz analyses of this problem (refer- 
ences 4 and 5). Furthermore, the orthogonality properties 
of the simple Fourier expansion lead to energy expressions 
of a simplicity that is instrumental in permitting accurate' 
computations. 

Approximate solutions: upper and lower limits. — 

The Lagrangian multiplier method, as well as the Rayleigh- 
Ritz method, gives a theoretically exact solution for the 
buckling stress; but ordinarily only approximate results are 
obtained because of the practical necessity of considering 
finite rather than infinite determinants. In the Rayleigh- 
Ritz method the approximate result is always higher than 
the true buckling stress. In the Lagrangian multiplier 
method, however, it is possible to obtain approximate solu- 
tions in two different ways, which permit the computation 
of a lower limit as well as an upper limit to the true buckling 
stress. As determinants of higher order are used to obtain 
approximations of higher order, both the upper-limit and 
lower-limit results approach the true buckling stress. Thus, 
the Lagrangian multiplier method can be used to provide a 
result to within any specified degree of accuracy. It may 
be expected, furthermore, that a particular determinant in 
the Lagrangian multiplier method ought to yield a more 
accurate result than a determinant of equal order in the 
Rayleigh-Ritz method. 


LAGRANGIAN MULTIPLIERS 

The procedure used in applying the fundamental mathe- 
matical principles of Lagrangian multipliers is described in 
this section; a general proof of the validity of the method is 
given in the appendix. 

Let it be required to minimize a function of N variables 

./Vl, X 2 , X;i, . . . X N ) (0 

where the .r’s are not independent but are bound together 
by the relationship 

*p(.r, , Xt, x 3 , . . . r v ) = 0 (2) 

Lagrange’s method of simultaneously minimizing / and 
satisfying the constraining relationship (2) is to minimize' 
the function 


J—\ip 

with respect to the .r’s. The quantity A is the undetermined 
Lagrangian multiplier. The necessary conditions for min- 
imizing/ then become 


- 0 

d.r K Ox K 


<p = 0 


(IT-= 1, 2, 3, 


AO 


(2) 


\ote that these expressions are A r +1 equations in the A + l 
unknowns x-,, . . . and A. 

If there' are two relationships that constrain the r’s; 
that is, if 

<P\ ( i , J '2 , - - • X \ ) 0 

<p-A- r i , L, • • • x N ) = 0 


two Lagrangian multipliers are then needed. The function 
to be minimize'd becomes 


f — A i<pi — A 2 V2 

and the minimizing equations are 


df 

0.r A - 


-Aj 


d.r A 


dtpi 

dx K 


(K= 1, 


. A/ 


Vi=0 

^2=0 

The method is easily extended to cover the case of any 
number of constraining relationships. 

PRELIMINARY ILLUSTRATIVE EXAMPLE 

Before the main example is given, a simpler buckling 
problem will be analyzed by the Lagrangian multiplier 
method in order that the method of application of Lagrangian 
multipliers may be most clearly presented without the ob- 
scuring details of analysis of more complicated problems. 
This elementary problem requires the use of but a single 
Lagrangian multiplier, which leads to a single stability 
equation. 
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Consider a square plate, clamped along two opposite 
edges, simply supported along the other two edges, which is 
loaded in compression on the simply supported edges. (See 
sketch.) From the exact solution of this problem, the deflec- 
tion surface of the plate is known to be sinusoidal in the 
x-direction. The deflection in the ^-direction, known to he 
symmetrical, must satisfy the clamped-edge boundary 
conditions; that is, zero deflection 


and zero slope 


w(r, 0) — w(r, a)— 0 


Oro 

dy 


(*< 0) (JT; ") <> 


(«) 

( 4 ) 


The present method uses a cosine-series expansion, whereas 
the Trefftz procedure would use a sine-series expansion. The 
problem is solved by both methods for comparison. 



Cosine-series solution. -In the cosine-series solution the 
expansion 


«>=sin 


m 7r x 

<’ Zj 

n= 0,2, I 


COS 


tlTTlJ 

a 


( 5 ) 


may represent, the deflection surface having m half-waves 
in the x-direction, since the Fourier series of even cosines is 
a complete symmetrical set. 

The boundary conditions (equation (4)) on the slope are 
satisfied by each term of the expansion; however, in order 
that w satisfy the conditions of equation (3) on the edge 
deflection, it is necessary that 



m=0,2,4 


Equation (6) is a constraining relationship on the a’s and 
as such will be introduced in Lagrange’s minimization process. 

As in the Rayleigh-Ritz method, the internal energy of 
deformation and the external work of the stresses are then 
calculated. Using the value for w as given by equation (5) 


in the general formulas (reference 6, equations (199) and 
(201) (modified)) 



gives 

l r -=g^ ^ (m' 2 + n l Y ( 1 + <5 n „) a 

11=0.2.4 

m i w 1 a r t \ ' . 

I— a / (l+5o«)a„" 

s Z_J 

11=0,2,4 


The usual 
expression 


Rayleigh-Ritz procedure 


V—T 


requires 


that the 

( 9 ) 


be a minimum with respect to the a’s. In the present ex- 
ample, however, the a’s are not independent but are bound 
by equation (6). Hence, mathematically stated, the ex- 
pression V- T must be a minimum, subject to the constraint 
relationship on the a’s 



11 = 0 , 2.4 


Solving this minimization problem by Lagrange’s method 
makes it necessary to minimize 

(F-D-X^o, (10) 

n =0,2,4 


with respect to the a’s. The necessary conditions for a 
minimum then become 



11 = 0 , 2,4 


or, upon differentiation and simplification, 

4 a 2 

(lFMK^+i 2 )*-™ 2 *]^- Jx) x =° (i=o,2,4, . . .) (ii) 



M-U.2,-1 
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Solving equation (11) for a } and substituting into equa- 
tion (6) gives the stability equation that determines k: 

[(m 2 H-J 2 ) 2 — m 2 k] (1+S 0 <)~ ° 

j= 0 , 2,4 

For a particular number of half waves m, this equation 
may be solved by evaluating the series for several trial 
values of k and interpolating to find the k that makes the 
series vanish. The correct value of m is that which gives 
the lowest value of k. For two half waves (m= 2) in the 
loaded direction, the theoretically exact value of £ = 7.69 
(reference 6, p. 345) is obtained when only 10 terms of 
equation (12) are computed. 

Sine-series solution (Trefftz method). — The same problem 
will be treated in the manner suggested by Reissner (refer- 
ence 2), which is similar to Trefftz’s method (reference 1). 

Let 

. mirx VV , . my , , ,,, 

a 'a v 7 

r=l,3, 5 

The boundary conditions on deflection (equation (3)) are 
now satisfied term by term, but the conditions on the edge 
slopes (equation (4)) are satisfied only by making 

V rb t = 0 (14) 

r=l, 3, 5 

Now, the expression V— T is computed from formulas (7) 
and (8) by using the value of w given by equation (13); 
then, by application of Lagrange’s procedure, 

{V-D-y^ rb T (15) 

r — 1 . 3, B 

must be minimized with respect to b„ where s=l, 3, 5, . . .. 
The Lagrangian multiplier is y. 

The minimization equations are 

4 a 2 

[ ( TO*+« i ) s -mm-^)7*= 0 (*=1,3,5, . . .) (16) 
^ rb T = 0 (14) 

r= 1,3, 5 

Solving equation (16) for b s and substituting in equa- 
tion (14) gives as the stability equation 

(m 2 -j-s 2 ) 2 — m 2 £ ^ 

*=1,3,5 

Comparison and discussion of results. — The series 
in equation (17) converges approximately as 1/s 2 , whereas 
the series in equation (12) converges approximately as 1/j*. 
Because of the more rapid convergence obtained in the 


stability equation, the Lagrangian multiplier method is 
preferably used to satisfy the zero-deflection condition rather 
than the zero-slope condition. Slope is the derivative of 
deflection, and, in general, differentiation of a Fourier series 
makes it more slowly convergent. 

THE COMPRESSIVE BUCKLING OF A RECTANGULAR PLATE 
CLAMPED ALONG ALL EDGES 

The previous elementary example required only a simple 
Fourier expansion and but one Lagrangian multiplier to sat- 
isfy the boundary conditions. The more difficult problem 
of finding the buckling stress of the rectangular plate 
clamped on all edges and loaded as shown in the accompany- 
ing sketch necessitates a double Fourier series, as well as an 
infinite set of Lagrangian multipliers to satisfy the boundary 
conditions. 



Boundary conditions. — -The boundary conditions of the 
problem are: 

Zero deflection, loaded edges 

w(0,y)=w(a,y)=0 (18) 

Zero deflection, unloaded edges 

w(x, 0)=w(r, b) = 0 (19) 

Zero slope, loaded edges 

g(0,V)=g(a,lf) = ° (20) 

Zero slope, unloaded edges 

§)!(*, 0) = | 1 y(*,b) = 0 (21) 

Fourier expansions. — In order to achieve a rapidly 

convergent solution, the principles established by the pre- 

ceding example are used as the basis for choosing the Fourier 
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expansions to satisfy, term by term, the conditions of zero 
slope rather than those of zero deflection. 

The buckling deformation corresponding to the lowest 
buckling stress is always symmetrical perpendicular to the 
direction of load but, depending on the aspect ratio of the 
plate, may be symmetrical or an Asymmetrical in the direc- 
tion of load. Thus, for symmetrical buckling, let 


Since this Fourier series must vanish, each infinite series that 
constitutes a coefficient of a cosine term must vanish. (All 
the Fourier coefficients of the Fourier expansion of the func- 
tion zero are zero.) Hence 

«*= 0 (3 = 0,2,4, . . .) (25) 

m= 0, 2, 4 


VA 


- X 

m= 0 , 2 , 

and, for an Asymmetrical buckling, let 


rnwx iiiry 
a mn cos a cos-j-'- 

m=0, 2, 4 n =0,2,4 


( 22 ) 


V V >'« n *y 
" 2j Zj "" COS a COS b W 

r=l, 3, 5 71=0.2,4 

It is sufficient, for purposes of demonstration, to consider 
only the case of symmetrical buckling. Hereinafter, w there- 
fore refers to the value given by equation (22). 

Energy expressions. — Using the expansion given by 
equation (22) in the evaluation of the general energy and 
work integrals of equations (7) and (8) gives 


w^Db 

8a s 




t rajb 
8 a 


X X (m 2 + n 2 fi 2 ) 2 { 1 + 8 m n + 8 0n )a„ 

11 = 0,2,4 n—0 f 2, 4 

X X /«.-( 1 + 8 0m ) a m n~ 


Wl=0,2,4 7? = 0,2,4 


Then 


where 


,T - 7 ’-s“ X X 

771=0, 2,4 77 = 0,2,4 


(24) 


1 

A„, 


= [ (m 2 + /A/3 2 ) 2 - km 2 ?} ( 1 + S m + 5 n „ ) 


Note that V— T is independent of a w , since 


Constraining relationships. — The boundary conditions 
of zero slope (equations (20) and (21)) are satisfied by 
each term of the expansion of equation (22), but the condi- 
tions on deflection (equations (18) and (19)) must be satisfied 
by the expansion as a whole. Substituting win to equation (18) 
gives, along the loaded edges, 


w(0, y)=w 


{a, y) — ^ rt,„ 0 + cos ^ ^ a m . 


771=0, 2, 4 


771 = 0, 2, 4 


+ cos 


4rr y 


X a - 


+-?X " 

)« = (), 2,4 

+ ■ ■ • =o 


By expressing the fact that there is zero deflection along 
the unloaded edges (equation (19)), it can be similarly shown 
that 

^ a in —Q (f=0, 2, 4, . . .) (26) 

n =0, 2, 4 


Now, U— T must be a minimum with respect to the a’s, 
which are bound by equations (25) and (26). As the prob- 
lem now stands, however, it is not in the form to which 
Lagrange’s minimization process can be applied since V—T 
does not contain a 00 , whereas the constraint relationships 
do contain « 0 o- Hence, u 0 o is eliminated from the con- 
straint relationships by subtracting the first of equations 
(26), the equation for i=0, from the first of equations (25), 
the ecpiation for j= 0. The final set of necessary con- 
straining relationships on the minimization of the energy 
(expression (24) 

03 

V -J_ 0 2 

/ A iL mn 
/ i -n-mn 
4 77=0, 2, 4 


V-T 


T [Db 

' 8a 3 


X 

7/7=0, 2, 


then becomes 


^ ^ ttmO ^ ^ & 0 t i — 
~ ‘ * 7 = 2, 4,6 

X “■ 

»=(), 2,4 

X “• 


777 = 2, 4,6 77 = 2,4,6 


= 0 0 = 2, 4, 6, ...) 


Y (27) 


0 (i=2, 4, 6, . . .) 


77 =0, 2.4 


Theory of upper and lower limit solutions.^! theoretically 
exact solution to the problem would be obtained if the 
energy expression (equation (24)) were minimized with 
respect to all the a’s and at the same time all the rela- 
tionships (27) were satisfied. This result follows from the 
facts that: (a) the expansion of equation (22) is a complete 
symmetrical set, capable of representing the exact symmet- 
rical buckling deformation, and (b) the fulfillment of the 
conditions of equations (27) ensures that the boundary con- 
ditions are completely satisfied. An exact solution is not 
possible, however, because it would involve an infinite deter- 
minant-, so that two different modifications of the ideal pro- 
cedure are used to obtain approximate results. One of these 
methods gives an upper limit to the true buckling stress, 
whereas the other gives a lower limit. 
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An upper limit to the buckling stress can be found by 
arbitrarily setting some a ’ s equal 'to zero, minimizing expres- 
sion (24) with respect to the remaining a’s, and satisfying all 
the constraint relationships (27). An upper limit is obtained 
inasmuch as arbitrarily setting some of the Fourier co- 
efficients equal to zero has tin' effect of restraining the deflec- 
tion of the plate, which in effect stiffens the interior of the 
plate and increases the stress required to buckle it. 

A lower limit to the buckling stress can be found by 
minimizing expression (24) with respect to all the a’s but 
satisfying only some of the constraining relationships (27). 
Neglecting some of the constraining relationships has the 
effect of giving the plate greater freedom at the edges and 
hence reducing the stress required to buckle the plate. 

Lower limit solution. — In accordance with the require- 
ments for a lower limit, the constraining relationships (27) 
will be satisfied only up to j—q and i=p. By Lagrange’s 
minimization process, the function to be minimized is then 


n __**Db V V 1 - / V 

0 ■ 8a :i 2 j Zj d,„ n " “ \ Aj a " 

■nt = 0, 2, 4 n = C), 2, 4 \m = 2, 4, (5 

U V 

- X a "‘ J ~ X X 

.#-2, 4, « m=i), 2, 4 ( = 2,4,1! n-O.2.4 



The a, X 2 . . . X s , H 2 . . . n,, are Lagrangian multipliers. 
The equations for minimizing V—T with the constraining 
relationships (27) on the a’s then become 


dG 

da„,„ 


0 , 2 , 1 , . 


Equations (27) up to i—p,j 



(28) 


Evaluation of dG/da,„„ gives 

^ -/ a mn- <x(8 0n —8 tn0 )—\ n —n m =() (29) 

where X„ appears only if 2^-n ^q and n m appears only if 
From equation (29) 


(Ipid , ft+e d (d, & ■ 2, 4, 6, . . .) 


Hence for any particular d, e, either 

1 

J d p\-d, 

or 

U pp <l . q i-e t) 


0 


The first alternative, however, ordinarily would require k to 
be very high, corresponding to the buckling stress of a 
buckling mode with many waves in both directions. For the 
lowest buckling load, then, 


<i +■« — d id, p 2, 4, G , . . .) 


It is therefore necessary to be concerned with oidy the other 
a’s, which, from equation (29), are 


4 ($ 

n : ^Y)[) M/«) 


4a 3 


^ mi) ( ^ ~i M/«) 


4a 3 . . 

0 n ~~ 4 Aim V ^ \ ^ n) 


(to, 0) (30) 


In equations (30), X„ does not. appear if //> q and n„, does 
not appear if to)> p. 


Substituting the values of a’s given bv equations (30) back into the constraining relationships (27) up to j=q, i=p gives 


A mo -f- ^ -In// J h * 1 i .... -lo//X„ — 0 

= 2,4,1! 11-2,4.8 / ill='2, 4, li ;i-2,4,(> 

P cc 

— A n] a + ^ A mjUm+Xj ^ A mj = 0 ( j— 2, 4 , 6 , . . . q) 

■m= 2,4,(» 7/1=0, 2, 4 

CO <1 

A n,a~\ i ^ d,»+ ^ A,„x„ = 0 ((=2,4,6, . . . p) 


(31) 


71 = 0,2,4 » = 2,4.1! 


These equations form a set of f(p + </) + 1 linear homo- 
geneous equations in a, fi 2 ■ ■ ■ t>„, X 2 . . . X 9 . iSince when 
buckling occurs the a’s are not all zero, by equations (30), 
the Lagrangian multipliers are not all zero. In order that 
equations (31) be compatible, the determinant of the coeffi- 
cients of the Lagrangian multipliers must vanish. The 
vanishing of this stability determinant provides the deter- 
minants! equation that may be solved for k by substitution 
of trial values and interpolation. 


That certain elements of the determinant consist of an 
infinite series of A mn terms is evident; these series converge 
rapidly. Since such rapidly convergent series are cal- 
culable. to any degree of accuracy, they may be considered 
as known quantities. Each value of A mn represents the 
potential-energy contribution of a term in the expansion for 
to; hence, the effects of infinite subsets of expansion terms 
enter into this solution. Thus, for p = q=2, the expansion 
terms corresponding to the a’s shown in table I enter into 
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the solution; similarly, the terms represented in table II 
enter into the solution when p — q — 4. 


TABLE I.— FOUR INFINITE STRIPS OF FOURIER 
COEFFICIENTS OF EXPANSION TERMS 


11 

0 

2 

4 

0 s ... 

ft 

fffto 


«40 

iiijii dsn 

2 

flf.2 

<1 22 

in i 

fl«2 «<*■> 

1 

lliM 

«24 



0 

I/O.', 

«2K 



s 

(his 





It is to be emphasized that the redundancy of a constrain- 
ing relationship is a peculiarity of only the upper-limit solu- 
tion, since, as shown by the proof given, the redundancy 
depends on the fact that all the constraining relationships 
must be satisfied. 

With the elimination of the redundant condition, the 
necessary constraint relationships become 

2 0 U ■■■-■■ '2, 4, «,...) (32) 

III — 0. 2, 4 


TABLE II.— SIX INFINITE STRIPS OF FOURIER 
COEFFICIENTS OF EXPANSION TERMS 




o 

-i 


ir.’ti 
< i-i- 


'to,- 


dpi 


a v 

dpi 


Hull 

a..< 


it,.i 


(ISO 


llu 



Upper-limit solution. -The lower-limit solution satisfied 
only some of the constraining relationships (27) but assumed 
the existence of all the Fourier coefficients. If an upper 
limit is to be obtained, it will be necessary to satisfy all the 
constraining relationships while arbitrarily assuming some 
a ’ s to be zero. 

As a direct result of the necessity of satisfying all the 
constraining relationships in the upper-limit solution, it is 
found that the first of equations (27) is redundant and may 
be discarded, since it is automatically satisfied when all the 
remaining equations (27) are satisfied. As a proof of this 
redundancy, the conditions 



7/i =0, 2, 4 


are summed over j and subtracted from the sum of the 
conditions 


^ <(i,u-<) (i* 2,4,0,...) (33) 

n = 0,2,4 

At this point, in accordance with upper-limit theory, it 
is necessary arbitrarily to set certain a’s equal to zero, i t is 
possible, to take advantage of the Lagrangian multiplier 
method by allowing infinite rather than finite sets of a’s 
to exist and still to obtain a stability determinant of finite 
order. Thus, infinite strips of coefficients of the type shown 
in tables 1 and 2 can enter into the solution. In the lower- 
limit. case, the existence of all coefficients was assumed, bid 
the coefficients a v+dwQ+e were proved to be zero; in this 
upper-limit solution, it will lie arbitrarily assumed that 
these same, a’s are zero; thus 

( ip+tt , 0+e — 0 ((h r 2, 4, h,... ) 

The constraining relationships (22) and (22) become 


V' t) 

JLJ 

m=o. >. ^ 


p 



L-! 

n—tt, 4 


( j -- 2, 4, 6, ...(/) 

O’ <1 F2, 7-f 4, ...«>) 
0' = 2, 4, 6, . . . p) 


(22a) 


(22b) 


(22a) 



W — 0 , 2 , -1 


over i and g*ive 

i = 2, 4,0 « =0, 2, 4 

Simplifying this equation 

®il) 

/ — 2, 4,0 




Tj t amJ '~ 


j = 2,4,6 m- 0,2,4 


o 

^7 <h,i~ 0 (/ =P + 2, p + 4, ...“>) (23b) 

n — 0, 2, 4 

The function to be minimized is 

a - y S 2 ) E X x - E - 

7ii =0, 2, 4 71 = 0,2,4 / = 2, 4, 6 m = 0, 2, 4 

00 P P ” 

^ ^ ^7 ^ ^ ^ ^ M? a in 

7 = i/+ 2, i/+4 714=0,2,4 f' = 2,4.0 71=0,2,4 

</ 

Me ^ 

i - /j+2, //+4 // -0, 2, 4 



which is precisely the first of equations (27). 


i 
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The first double summation of this equation extends over 
only the values of m and n such that 
mS p if +></ 

'R = 2 if m)>p 

Setting < — = 0 then gives, for all the a’s arbitrarily allowed 

00 mn 

to exist, 

A.o? 

dm u — 

where X 0 and ju t) do not exist. Substituting back into the 
constraint equations (32a), (32b), (33a), and (33b) gives 


x ' 2 
m= 0, 2, ■ 

4 .+ 
* 1 

4 

£ ,w.,= 

m = 2, 4, 6 

0 (i= 2,4,6, . . 

■ q) 

(34a.) 

P 

?/i=0, 2, - 

A mi -\- 

1 

P 

A m j fx fa — 

7/1 = 2, 4, 0 

0 {j=g_+2,q+4,. 

■ ■ co ) 

(34b) 

t 

A OiXn + Mi 

^ A in = 0 

(7 = 2, 4, 0, . . . 

■P) 

(35a) 

« = 2, 4, f 


71=0, 2, 4 




( l 

2 

A. i n \ n ~ ix i 

<7 

^ ^ A fn— 0 

ti=p + 2,p+4, . . 

■ 03 ) 

(35b) 


fl = 2,4,6 71=0,2,4 


These equations involve all the Lagrangian multipliers. 
They can be reduced to a set of equations, however, in 
X 2 . . . X„ g 2 ... up in the following manner: 

From equation (35b), for i=p-\- 2 , p + 4, . . . 

<? 

. 1 i n X a 

71 = 2,4,6 

A tn 

71=0,2,4 

From equation (34b), for j=q-\-2, 2 + 4, . . . 


p 



7/1=0, 2, 4 


Substituting these expressions for /u* and X, into equations 
(34a) and (35a), respectively, gives as the final stability 
equations: 




+ J A Mj + J 


A m i 


r ^ 

n = 2,4,6 


?/l=0, 2,4 m- 2,4,0 


m=p- 1-2, />- M 


<7 **2 

2 

= 2, 4, € 

Q 




f/ CO 

^ ^ ^ ^ A- nl 


n = 2, 4, 0 


7/ =<7+2, </-f4 


S- 

m=2, 4, 6 
P 

S 


mn 

V 11 = 072,4 > 


= 0 (i=2,4,fi, 


mnP-m 


A n 


m n 

^ m =072, 4 


+m< Aj„=o 

71=0,2,4 


ff) 


0'=2, 4, 6, 


7>) 


(36) 


Equations (36) form a set of ~ (p — //) linear homogeneous 

equations in X 2 . . . X,„ /i 2 . . . /i p . The stability determinant 
is the determinant of the coefficients of the X’s and +S. 

It is of interest to note that in the usual Rayleigh-Ritz 
solutions only finite sets of expansion terms are ever taken 
into account, and the order of the determinant obtained is 
ordinarily equal to the number of terms considered. It is 
then reasonable that a particular determinant obtained by 
the Lagrangian multiplier method, which considers infinitely 
more expansion terms than a Rayleigh-Ritz determinant of 
equal order, may be expected to give a more accurate result. 

Numerical example. For the case of a square plate, 0=1, 
upper and lower limits were computed. The results for- 
th e buckling-stress coefficient k were: 


Lower limit j Upper limit j 
, 

9.99 1(1.11 

10.07 ; 10. OH j 


The expectation that the Lagrangian multiplier method 
should give closer upper limits than the Rayleigh-Ritz 
method, for a given-order determinant, can be confirmed for 
this example. A second-order Lagrangian multiplier deter- 
minant gives an upper limit of&=10.11, whereas Maulbetsch 
(reference 4) and Smith (reference 5) use complicated de- 
flection functions in the Rayleigh-Ritz method to derive 
third-order determinants that give, respectively, k= 10.48 
and k= 10.11. 

It is seen that the second approximation, requiring the 
evaluation of a fourth-order determinant for the upper limit 
and a fifth-order determinant for the lower limit, definitely 
establishes the value of k to within 0.1 percent. 

Levy (reference 3) used an ingenious method of obtaining 
lower limits that is, in fact, equivalent to the Trefftz method 
of using double sine series and satisfying the zero edge-slope 
condition by the Lagrangian multiplier method. On the 
basis of computations involving determinants up to order 
twenty, Levy concluded that his results obtained from tenth - 
order determinants are within 0.1 percent of the true results. 


Approximation 


First; j>=//=2 
Second; p=q= A 
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Inasmuch as Levy obtained 4=10.074 for the square plate, 
the present relatively simple upper- and lower-limit calcula- 
tions show that his estimated limit of error is correct for 
this case. 

CONCLUSIONS 

1 . The Lagrangian multiplier method can be used to 
compute accurate upper and lower limits to the compressive 
buckling stress of a clamped rectangular plate, thereby 
bracketing the true result. 

2. From a consideration of rapidity of convergence toward 
the exact solution in clamped-plate problems, it is preferable 
to use an expansion that satisfies the zero-slope boundary 
conditions term by term rather than the zero-deflection 
boundary conditions. 

3. Because of the fact that the Lagrangian multiplier 
method permits the effects of infinite subsets of expansion 


terms to enter into the solution, it is believed that a particular 
stability determinant derived by the Lagrangian multiplier 
method will, in general, yield a closer upper limit than that 
obtained from a determinant of equal order in the Rayleigh- 
Ritz method. 

4. It is expected that the method of Lagrangian multipliers 
may be useful in the analysis of other stability and vibration 
problems. In particular, the method may be immediately 
applied to the determination of vibration frequencies of 
clamped plates, and to the determination of buckling stresses 
of clamped plates under compression in two directions. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., May 3, 1946. 



APPENDIX 


GENERAL PROOF OF THE METHOD OF LAGRANGIAN MULTIPLIERS 


Let it- be required to minimize 

/(X 1, X 2 , X ? „ . . . X v ) 


(Al) 


where the N x ’ s are bound by the P independent relation- 
ships (P<fN) 

<Pj (x u x 2 , x-i, . . . 3w) = 0 (J= 1, 2, 3, . . . P) (A2) 
Tt will be proved that the equations for determining the 


minimizing values of the x’s are: 


df , dp 


dpi 


dp 


dx, 


l+X^r+X^+.-.+X,^ 0 {K=\, 2,3,... A) (A3) 


' dx 


dx, 


<pj {x u x 2 , x 3 , . . . x. v ) — 0 (J= 1, 2, 3, . . . P) 


(AS 


The X’s are Lagrangian multipliers; these A T-\-P equations 
determine N x’s and P X’s. 


If the values of only N—P x’s are known, the remaining P x’s are determined from the P relationships (A2). For 
convenience, consider the last P x’s in equation (Al) to be dependent upon the first N—P x’s. Then for f to be a 
minimum its first partial derivatives with respect to the independent x’s must vanish, or: 

df dxv_pf2 


df j df dx x . r . : . 

dx M dx N _ F+1 dx M dx A r_p +2 


df dx. v ^ 
dx M 1 dx A - dx M 


(M= 1,2,3, . . . (N-P)) 


(A4) 


But each of these equations contains P quantities that cannot be directly evaluated -the derivatives of the depend- 
ent variables with respect to the independent variables. For each value of M, these P derivatives are determined by 
differentiating each of the P constraint relationships (A2) with respect to x M . Thus, 

dXy_ P+ 2 


dtpj 


dyy, di pj dx.v- /-.a i _j_ 

dx M dx, v /■(. i dx M dxy^p^ a 


dx M 




, dtpj dx.v „ 

^ dXv dX;, 


b/ - 1,2,3, . . . P) 


(AS) 


Now, for each particular value of M, equation (A4) and the P equations (Ao) make up a set of P-fl linear homo- 

Sinco these quantifies are surely not all zero, 


, • • i, i , , , dx v _p_i dxjv_p +2 dx. v 

geneous equations m the P + 1 quantities 1, v • ■ > ■ ■> . . . < 

OX \f OX\f OX nf 


be true that: 


ients 

must v 

anish. 

Hence, 

I dj 

d / 

d f 

df 

dx„ 

dx.v_p+i 

C)X\t_ P42 

dx,\ 

dp, 

dy, 

dv?i 

dip, 

d.r„ 

d.x.v p 1 1 

dXv ..p-\ 2 

dx A 

d p-> 

d p-> 

d p-> 

dp- 

dx AI 

dx.v-. p+i 

dX.v-p.f2 

dx.x 

d pp 

d <pp 

dip,- 

dp, 

dx*, 

dx.v_pj.i 

d/y- /■» -j-O 

dx.y 


(A6) 


It will now be demonstrated that these necessary minimization equations will hold if equations (A3) hold. Interchanging 
the rows and columns of the determinant in equation (A6) gives: 


df 

dp. 

dp 2 

dp,- 

dx M 

dx.M 

dx.„ 

d.r.u 

df 

dp, 

dv32 

dp,- 

dXy_ p ,_i 

dx.v.. p. |-i 

dx.v-p.,, 

dx,v-_p.j. i 

df 

dy>i 

dp-> 

d^p 

dx.v p 

dx.v_p.f2 

C)X^y_p j_2 

dX.v-p + 2 

df 

dp. 

dpi 

dpp 

dx N 

dx.v 

dx.v 

dx.v 
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Tlu* vanishing of this determinant is, however, precisely the condition of compatibility of the equations 


d f 
d.r„ 

, x d^i 
^ X| d.r v 

+ X 2 d.rJ +• 

i x £><p7 ’ 
'+ X "0x,, 

> 

-- 0 

d f 

d.r,v r f i 

+ X| d,w. 

_L x d ^-' 4- 

+ -dr + ‘ 

V' l N -r h 

dXtf-p 4 1 

= 0 

Of 

d.r. v n 


+ x -' d , d<p ' +■ 
U - , A ! -P\2 

• d A,., 5 ^ 
d.r.v -P 4 2 

— 0 - 

Of 

Ox v 

+ 

, x d{p 2 4- 

r X " d.r.v 1 ' 

.+ 

ipX N 

= 0 


when they are considered as linear homogeneous equations in the quantities 1, Xj, X 2 , . . . \ P - 

Since a determinant (A7) exists for each value of M up to (N-P), a set of equations (A8) exists for each M. It is 
seen that in these sets only the first equation varies, since only the first equation depends upon M. Observation 
shows that all the (N-P) determinants of equation (A7) can be derived from the set of A 7 equations 


df 

0x K 


+ Xi 


dg>i , > 
d.r A r ^0x K 


+ 


■ + X,> 


£><p/ 

dr a 


(K= 1,2,3, 


AO 


(Ad) 


by successively writing the determinants of compatibility 
of the last P equations with each of the first (N—P) equa- 
tions in turn. It has thus been proven that if equations (Ad) 
are true, the minimizing equations (AO) , equivalent to 
equations (A7), must hold. 

It is seen, however, that equations (Ad) are N equations 
in (N+P) unknowns consisting of A 7 x’s and P X’s. The 
remaining necessary P equations come from the original 
equations of constraint (A2). Hence, the simultaneous 
equations (A2) and (Ad) 


d/ i , d<pi , dv ? 2 , ,, 

d/ A + Xl dr A +X 2 dx A + •■• +Xj 


d tpp 

0x K 


= 0 (K= 1,2,3, 


• N) 
(A3) 


<pj(A, A>, r 3 , . . . r v i=0 (J= 1, 2, 3, . . . P) (A2) 


are necessary equations for the minimization of /(*,, x- 2 , x :s , 
. . . x N ), which was to be proved. 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Longitudinal . 

X 

X 

Rolling __ . 

L 

Y >Z 

Roll 


u 

V 

Lateral 

Y 

Y 


M 

Z *X 


e 


q 

Normal 

Z 

Z 

Yawing _ 

N 

X >F 

Yaw - 


10 

r 











Absolute coefficients of moment 

p ft M p . AT 

Vl ~qbS Cm ~qcS ° n ~qbS 

(rolling) (pitching) (yawing) 


Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D 

V 

Diameter 
Geometric pitch 

P 

p 

Power, absolute coefficient C P — — 

’ pn 1 * 3 D B 

P/D 

V' 

Pitch ratio 
Inflow velocity 

a 

Speed-power coefficient 

v. 

Slipstream velocity 

V 

Efficiency 

T 

Q 

• T 

Thrust, absolute coefficient C T = — 

pjtSLr 

Torque, absolute coefficient C Q — tfjja 

n 

4> 

Revolutions per second, rps 
Effective helix angle = tan~ 


1 hp =76.04 kg-m/s=550 ft-Ib/sec 
1 metric horsepower =0.9863 hp 
1 mph= 0.4470 mps 
1 mps=2.2369 mph 


5. NUMERICAL RELATIONS 

1 lb =0.4536 kg 

1 kg=2.2046 lb 
1 mi= 1,609.35 m=5,280 ft 

1 m=3.2808 ft 



